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Abstract 

We study expectation values related to the dynamical lattice spacing that occurs in 
the recent 2d lattice super- Yang- Mills constructions of Cohen et al. |hep- lat /0307012 . 
The corresponding observable in the fully-quenched ensemble would appear to indicate 
a difficulty with the proposed continuum limit. However, we find that the same ob- 
servable in the phase-quenched ensemble takes a very different, perhaps encouraging, 
average value. Unfortunately, we are not able to obtain results for the full theory, due 
to the nearly flat distribution of the complex phase of the fermion determinant in the 
phase-quenched ensemble. 
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1 Introduction 



In [I], a lattice action has been proposed by Cohen, Kaplan, Katz and Unsal (CKKU) for 
the (4,4) 2d U(k) super- Yang-Mills. 1 The Euclidean target theory action is 



where s M (// = 0,1,2,3) are hermitian scalars, F is the 2d Yang-Mills field strength, and 
ipi (i = 1,2) are 2d Dirac fermions, all in the adjoint representation of U(k). In |2j, some 
aspects of the fermion determinant were examined, with results very similar to those in 
We will have occasion to review those results below. 

Here we will study a rather fundamental aspect of the construction: the proposed emer- 
gence of a "dynamical" lattice spacing. CKKU make use of an N x N lattice that contains 
link and site fields that are k x k matrices. The classical lattice action contains many 
zero-action configurations. CKKU expand the classical lattice action about a particular 
class of zero-action configurations (discussed below) that are characterized by a parameter 
a. For this reason we will refer to such a configuration as an a- configuration. In the limit 
a — > 0, N — > oo, the classical action tends to 1)1.1)1 . Thus a is interpreted as a lattice spacing. 
However, it is dynamical as it has to do with a particular background configuration for the 
lattice fields. This strategy is based on the ideas of deconstruction jU E|. Studies where 
a partial latticization of 4d supersymmetric theories has been obtained by this approach 
include [HIIZI- 

The validity of the semi-classical expansion, about an a-configuration, rests on the as- 
sumption that it gives a good approximation to the behavior of the full lattice theory. But 
all a-configurations are energetically equivalent. Furthermore, there exist other zero action 
configurations that do not fall into the a-configuration class (shown below). Finally, it is 
important that the fluctuations about the a-configuration be in some sense subdominant. 
For these reasons, CKKU deform the action by adding an a-dependent potential that favors 
the a-configuration. The "continuum limit" then includes sending a — > in this potential. 
Although this deformation breaks the exact lattice supersymmetry, it is rendered harmless 
by scaling the relative strength of the deformation potential to zero in the thermodynamic 
limit. For this reason it has been argued by CKKU that the quantum continuum limit is 
nothing but the target theory. Based on the symmetries of the undeformed theory, it would 
appear that the target theory is obtained without the need for fine-tuning. For further 
details, we refer the reader to [T], as well as the articles leading up to it (HUH]- 

The deconstruction method for latticizing a continuum target theory does not require 
fermions or supersymmetry. The bosonic system (i.e., setting all lattice fermions to zero) 

1 The (4,4) 2d U(k) super- Yang-Mills is best defined as the dimensional reduction of AT — 1 6d U(k) 
super- Yang-Mills; the notation "(4,4)" denotes the number of left and right 2d chirality supercharges. 




(1.1) 
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contained in the model of CKKU already has the interesting feature of a deconstructed lattice 
Yang-Mills. The same semi-classical arguments yield U(k) Yang-Mills with 4 adjoint scalars; 
that is, the bosonic part of (jl.ljl . In the present article we will investigate the validity of the 
semi-classical argument in both this bosonic theory and in the full supersymmetric theory. 
We will take into account quantum effects, estimating key expectation values by means of 
Monte Carlo simulation. 

We now summarize the content and results of the rest of the article: 

• In the next section we introduce the essential features of the bosonic part of the CKKU 
construction that will be needed for the subsequent discussion. We then study the 
classical minima of the undeformed and deformed actions in Section El We show that 
in the undeformed theory there is a vast number of zero-action configurations. We 
further show that the deformation leaves only a particular a-configuration with zero 
action, modulo gauge equivalences. 

• In Section0]we outline the methods and results of our lattice simulation for the bosonic 
system. We find that the expectation values do not agree with the semi-classical 
predictions based on energetic arguments. In Section E] we suggest an interpretation of 
these results in terms of the analysis of Section El We find that the simulation results 
are consistent with entropic effects dominating the expectation values. 

• In Section El we attempt to include the effects of fermions. We show that this is 
crucial for understanding the dynamical lattice spacing in the supersymmetric theory, 
because the fermion matrix is highly singular for the zero-action configurations of the 
undeformed theory. Difficulties are encountered because of the complex phase of the 
fermion determinant. We present results for expectation values in the phase-quenched 
ensemble. However, we are unable to obtain reliable expectation values for the full 
theory, via phase reweighting, due to the nearly flat distribution of the complex phase. 

• In Section [7| we make some concluding remarks. In the Appendix, we provide a brief 
review of the M = 4 4d super- Yang-Mills moduli space, which arises in the discussion 
of Section El 

2 Quiver lattice Yang-Mills 

For CKKU, the starting point is the Euclidean M = 1 6d U(kN 2 ) super- Yang-Mills. In 
our considerations of the bosonic theory, the starting point will be just the Euclidean 6d 
U(kN 2 ) Yang-Mills. In either case, the action is dimensionally reduced to Od to obtain a 
U(kN 2 ) matrix model. The matrix model naturally possesses SO (6) Euclidean invariance. 
Next we note 50(6) D SO (2) x SO (2) D Z N x Z N . CKKU have identified a homomorphic 
embedding of Z N x Z N into the U{kN 2 ) gauge symmetry group. Using this, a Z N x Z N 
orbifold projection is performed to obtain a U(k) N Od quiver, or, product group theory. 2 In 

2 Quiver theories were originally studied many years ago in other contexts [jT)HTTj . 
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every respect we follow CKKU, except that in the bosonic theory we have set all fermions to 
zero. For details of the matrix model and orbifold projection, we refer the reader to pQ. In 
the interests of brevity, we will only give the final results. The undeformed bosonic lattice 
action is 
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The fermionic lattice action will be given in Section El where we will examine the effects 
of the fermion determinant. In 1)2.1)1 . x m ,y m ,z m are bosonic lattice fields that are k x k 
unconstrained complex matrices; m = (777.1,777.2) labels points on an iV x Af lattice, and 
i = (1,0), J = (0, 1) are unit vectors. The U(k) N symmetry is nothing but the local U(k) 
symmetry of the lattice action So, with link bosons x m in the i direction, link bosons y m in 
the j direction, and sites bosons z m , all transforming in the usual manner: 

Canonical mass dimension 1 is assigned to x m , y m , z m , whereas g has mass dimension 2. 

Although 5*o is a lattice action that describes a statistical system with interesting features, 
it is not in any obvious way related to a 2d continuum field theory (jl.l)) . As will be explained 
below, Sq > and a vast number of nontrivial solutions to So = exist, not all of which are 
gauge equivalent. In fact, the space of minimum action configurations, or moduli space, is 
a mult i- dimensional noncompact manifold with various branches (classes of configurations). 
For this reason it is difficult to say what a "continuum limit" might be; for there exists an 
infinite number of energetically equivalent configurations about which to expand, not all of 
which are gauge equivalent. 

A surprising result — pointed out by CKKU, and based on ideas from deconstruction — is 
obtained if one expands about the a- configuration 

x m = — ^1, y m = — ^1, z m = 0, Vm (2.3) 
ay2 ay2 

keeping g 2 = ga and L = Na fixed, treating a as small. (It is easy to see that So = for 
this configuration.) That is, we associate a with a lattice spacing (mass dimensions -1), even 
though it arises originally from a specific background field configuration. In this case, one 
finds that the classical continuum limit is nothing but the bosonic part of (Jl. 1)1 . which is a 
variety of 2d U(k) Yang-Mills with adjoint scalars. In the case of the supersymmetric quiver 
theory of CKKU, where fermions are present, one obtains (jl.l)) in full; i.e., (4,4) 2d U(k) 
super- Yang- Mills. 
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The trick is how to make the configuration ()2.3|) energetically preferred without destroy- 
ing all of the pleasing symmetry properties of the theory (This is particularly true in the 
super symmetric case.) In the quantum analysis, we must address the more delicate compli- 
cation of entropy as well. CKKU suggest a deformation of the bosonic action in an effort to 
stabilize the theory near the a-configuration ()2.3|) : 



S B 
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(2.4) 
(2.5) 



Here the strength of the deformation is determined by the quantity /i, which has mass di- 
mension 1. It is clear that the configuration (J2.3)) minimizes Ssb- (Other configurations that 
minimize So and Ssb will be discussed below.) Unfortunately, in the supersymmetric ver- 
sion of CKKU, the deformation S'sB breaks the exact supersymmetry of their original lattice 
action (hence the subscript "SB" = Symmetry Breaking). For this reason they demand that 
the strength of Ssb relative to So, conveyed by fi 2 , be scaled to zero in the thermodynamic 
limit. Thus we are interested in the effects of the deformation subject to this scaling. 

In much of what follows we will specialize to the case of U(2). This is merely because it 
is the simplest case and the most efficient to simulate. In this special case, x m ,y m , z m will 
be unconstrained 2x2 complex matrices. 



3 The classical analysis 

Here details are given of the classical analysis of the minima of the action Sb- In Section 
13. II we consider the undeformed action So', then in Section l3~2*l we consider the modifications 
induced by the deformation Ssb, which has the effect of lifting some flat directions in moduli 
space. Through understanding this classical picture, naive expections of what will occur in 
the quantum theory, based on energetics, can be formulated. 



3.1 Undeformed theory 

Here we neglect Ssb and examine the minima of So- Note that (|2.1|) is a sum of terms of 
the form Tr AA> (the first line involves squares of hermitian matrices). Thus So > with 
5*0 = iff the following equations hold true: 



XL ;X 



•^nl/n+i VnXn+j Vn^n+j ZnVn ^n-^n -^n^n+i (3-2) 

together with the h.c. of ()3.2|) . The set of solutions is the moduli space of the undeformed 
theory. 
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3.1.1 Zeromode branch 



To begin a study of the moduli space, we isolate the zero momentum modes: x n = x Vn, 
etc. Then Eqs. (|3.1|) and (|3.2|) reduce to 

[x\x\ + \y\y} + [z\z] = 

[z,y] = \y,z] = [z,x] = o (3.3) 

together with the h.c. of the second line of ()3.3|) . Eqs. ()3.3|) may be recognized as nothing 
but the D-flatness and F-flatness constraints that describe the moduli space associated with 
the classical scalar vacuum of M = 4 4d super-Yang-Mills. 3 The equations are invariant with 
respect to the global gauge transformation 

x — > axoc\ y — > aya\ z — >■ aza^ (3.4) 

Then it is well-known that solutions to ()3.3j) consist of x, y, z that lie in a Cartan subalgebra 
of U(k); the proof is reviewed in Appendix A. The global gauge transformations (|3.4jl allow 
one to change to a basis where this Cartan subalgebra has a diagonal realization. Thus one 
can think of the moduli space as the set of all possible diagonal matrices x, y, z, and all 
global gauge transformations ()3.4j) of this set. 

In particular, the zeromode moduli space of the undeformed U(2) theory is completely 
described by 

x = x° + x 3 a 3 , y = y° + y 3 a 3 , z = z° + z 3 a 3 , (3.5) 

with arbitrary complex numbers x°, x 3 , y°, y 3 , z°, z 3 , together with U(2) transformations of 
these solutions. 

Eqs. fl3.1j) and ()3.2j) also have non-zeromode solutions. We do not attempt to present 
an exhaustive account of them. We will merely point out a few such branches in order 
to illustrate that the undeformed theory has a very complicated and large set of So = 
configurations. This observation will be relevant to our interpretation of the simulation 
results in Section 

3.1.2 x m = y m = non-zeromode branch 

We have the very "large" branch of moduli space described by 

x m = y m = 0, z m = z Q m + z^a 3 , Vm (3.6) 

Again, zj^, are arbitrary complex numbers. Furthermore, site variable and thus 

transforms independently at each site as 

z m > Oi m z m ot^ m (3.T) 

3 We thank Erich Poppitz for pointing this out to us, as well as the branch of moduli space H3.6|l given 
below. 



5 



It can be seen that this branch affords a vast number of solutions to (j3.1j) and ()3.2|) ; 
there are N 2 such solutions, modulo choices for z^, G C and gauge equivalences. We 
will argue below that an understanding of the entropic effects that result from this branch is 
necessary to understanding expectation values in the quantum theory — even when a potential 
is introduced which gives these configurations nonvanishing action. 



3.1.3 z n = non-zeromode branch 

Another branch in moduli space is the following. First we set z n = 0, Vn, and introduce 
Fourier space variables 

Xn = a7 E ^'"^ V" = ^ E ^ k ' n ^> w = exp(27ri/iV) (3.8) 

k k 

where k = (fci,/^) an d &i>&2 £ [0,1, . . . , N — 1]. Then taking into account z n = 0, the 
conditions ()3.1|) and ()3.2|) are equivalent to: 

k 

= ^(w^Vw-^Vw/k) (3-9) 

k 

for all £ = (£ 1; £ 2 ) an d £1, £ 2 £ [0, 1, • • • , N — 1]. Next we turn off all modes except one for 
both / k and g k : 

/k = 4,k'/k', = 4 ,-k'P-k' (3.10) 

Here and below, no sum over k' is implied. When substituted into ()3.9|) . only 2 nontrivial 
conditions survive: 

= [fl, f v ] + [gl k „ g-u], = f v g^ - ^ +3 >* g- V fv> (3.11) 

For the U(2) case, we find that solutions exist if u>^ + ^' k = ±1. We already know from 
the zeromode considerations that for 

u (t+frV = 1 we have solutions for /k',<7-k' diagonal 
matrices. In the case of a/ z+ ^" k = —1 it is easy to see that there are solutions, say, of the 
form 

/k' = z f a 3 , g^ = z^a 1 + ba 2 ), z f , z g e C, b e R (3.12) 
There are many values of k' for which u/ z+ ^' k ' = ±1. For N even these are 

N 3N 

K + k' 2 = 0,-,N,— (3.13) 
For N odd, k[ + k' 2 = 0, N are allowed and in the cases where 

K + ^H±l,^pl ,3,4) 
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(J3.12JI yield approximate solutions to (|3.11|) . with an error of order 1/N. Thus in the N oo 
limit the number of Sq = configurations in this class is vast; in fact, it is easy to check 
that the number of such configurations is approximately 2N, modulo gauge equivalences and 
various choices for the constants in (j3.12|) . 



3.2 Deformed theory 

Now we consider the supersymmetry breaking deformation Ssb introduced by CKKU. To see 
its effect it is handy to rewrite the quantities that appear in it. Recall that complex 
2x2 matrix. Dropping the subscript, we can always define 



x = x° + x a a a , 



x f = x° + x a a a 



(3.15) 



Then it is straightforward to work out (// = 0, . . . , 3) 

xx ] = x^ + (x ^ + x°x c + ix a x*e abc )a c 



(3.16) 



Note that (j) x ^ are real, and that <j) x ' is positive definite. With similar definitions for 
the CKKU deformation is 
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(3.17) 



It can be seen that the deformation drives 0^ a , </>j^ a , toward the origin, and 0^°, 
toward l/2a 2 . When = 0, it is easy to see that = identically. 

To continue the analysis, it is convenient to rescale to dimensionless quantities using the 
parameter a: 



g = ga , /j, = /id, 



.2 ix 



ax r 



etc. 



(3.18) 



Then 
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(3.19) 



For any value of the lattice spacing a, the minimum of Ssb ls obtained iff 
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(3.20) 
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The conditions involving; just 

^m X m ~ 2' %m X m ~l~ X m%m ^m X m^ ~ ^ (3-21) 

Let us examine what additional constraint this places on classical solutions to S = 0, beyond 
the restrictions of the undeformed theory. 

First we note that neither of the non-zeromode branches discussed in Section 13.11 above 
are minima of Ssb- Thus we pass on to the zeromode configurations ()3.5|) . Eqs. 1)3.21)) then 
imply that (|3.5|) is restricted to the form 



x 



V2 



diag(e^,e-^) (3.22) 



and global gauge transformations of this. That is, x is restricted to be an element of the 
maximal abelian subgroup U(l) 2 of U(2), up to an overall factor of 1/V2. Similarly, we have 
for y, 

£ = _diag (e^,e-^) (3.23) 

Finally, 0^° = implies -2^m = 0> which has the unique solution z m = 0, Vm. 

Apart from global obstructions that are essentially Polyakov loops in the i or j directions, 
the configuration (|3.22|) and ()3.23|) can be gauged away. It is straightforward to verify that 
the required gauge transformation is (|2.2j) with 

a mi ,m 2 = e i(mi7 * +m27y) x diag [ e ^ m ^+ m ^y) , e -«(mi^+m 2 ^)^ ^ ^ 

This sets all x m , y m to unity except at the "boundaries" : 



£ iN-y x 



x N , m2 = __diag (e liV ^,e-^), Vm 
V2 



2 



e iNj v 



y mi ,N = ^-diag (e l7V ^,e-^), \/m x (3.25) 

For most purposes, we do not expect such vacua to distinguish themselves from the trivial 
vacua in the thermodynamic limit. In any case, global features such as these are typical of 
classical vacua of other lattice Yang-Mills formulations, such as the Wilson action. In our 
simulation study we will avoid this issue by restricting our attention to the expectation value 
of a quantity that is independent of these angles. 



4 Simulation of the bosonic theory 

The emergence of the effective lattice theory relies upon the assumption that fluctuations 
about this classical minimum are small, and that the equations ()3.20j) are a good approxi- 
mation to the corresponding expectation values in the quantum theory. For this reason we 
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study 



(e°} = (^J = (^Tr(x m xL)) (4.1) 
in our simulations, and compare the expectation values to the classical prediction (|3.20j) . 

4.1 Scaling 

We study (j4.1|) along a naive scaling trajectory: 

g 2 = a~ l g(a) = fixed (4.2) 

That is, we hold the bare coupling in physical units, #2, fixed; this is equivalent to neglecting 
its anomalous dimension. The dimensionless bare coupling g is then a function of a that 
vanishes linearly with a as the UV cutoff is removed. 

With regard to \x we follow the instructions of CKKU: we send the dimensionless coeffi- 
cient p, of the deformation Ssb to zero as 1/N while increasing N. 

pT 1 = cN, c = 0(1) (4.3) 

This is equivalent to scaling \i = 1/cL, where c is a constant and L = Na is the extent of 
the system. 

In the rescaled variables (|3.18jh the coefficient of the undeformed action is 1/g 2 , whereas 
the coefficient of the deformation is fi 2 /g 2 , as can be seen from (|3.19j) . Thus it is that 
the relative strength of the deformation vanishes in the thermodynamic limit, when ()4.3|) is 
imposed. 

We perform these scalings for a sequence of decreasing values of a. That is, we study the 
thermodynamic limit for fixed values of a. We then extrapolate toward a = to obtain the 
continuum limit. 

The physical length scales are set by g% and the system size L = Na. To keep dis- 
cretization effects to a minimum we would like to take g^ 1 ^> a. Equivalently, g^ 1 ^> 1. 
On the other hand we are most interested in what happens at large or infinite volume. To 
render finite volume effects negligible would require g^ 1 L. Equivalently, g^ 1 N. In 
the simulations we study the system for various choices of g^ 1 and N. We extrapolate to 
the regime 1 g~ l N, but often violate the bounds 1 < g^ 1 < N for specific points 
where measurements are taken. The reason for this is that data outside the optimal window 
1 g^ 1 «C iV is informative to the extrapolation. 

4.2 Sampling procedure 

We update the system using a multi-hit Metropolis algorithm. We attempt to update a 
single site or link field 10 times before moving to the next, with an acceptance rate of 
approximately 50 percent for a single hit. We find that this minimizes autocorrelations while 
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maintaining program efficiency. We have examined autocorrelations and the dependence of 
our observables on initial conditions. These studies have led us to make 500 themalization 
sweeps after a random initialization, and 100 updating sweeps between each sample. 1/2 to 
2 percent standard errors result from accumulating 1000 samples at each data point. 

4.3 Results 

In Fig. Owe show (0^°) [cf. (14.1)) ] as a function of N for various values of g~ l , having set 
c = 1 in (j4.3j) . Doubling g~ x is equated with halving the lattice spacing, according to ()4.2|1 . 
It can be seen that (0^°) tends toward smaller values as the lattice regulator is removed, 
contrary to the classical expectations ()3.20|) . At large enough values of N the curves flatten 
out to a constant value. 

To understand this behavior, we first note that we are computing an expectation value 
that is already nonvanishing in the undeformed theory (/i = 0). The undeformed theory 
expectation values 

OrW) = (r 1 ^ ) (4.4) 

for various values of N are shown in Fig. El The rescaling by g" 1 = a?g~ x is useful because 
it corresponds to removing g from the undeformed action Sq by a rescaling of the lattice 
variables [cf. (|2.1|) ]; this amounts to studying the undeformed theory in units of ^fg. It 
can be seen from Fig. 121 that f)4.4j) is far from zero, and is rather insensitive to N. The 
undeformed theory only contains two length scales, 1/yfg and N/^fg. For large N, it is 
not surprising that the local expectation value ()4.4|) is insensitive to the long distance scale 
N/ y/g. Rather, it is determined by the short distance scale 1 / J~g. The deformed theory 
retains this short distance scale. Once the system size is much larger than this scale, the 
local expectation value becomes independent of the system volume; this is particularly true 
because the relative strength of the deformation is being scaled to zero [cf. 1)4.3)1 ] . 

In Tabled we show the large N expectation values (|4.4|) in the deformed theory as well 
as those of the undeformed theory. Here, the values for iV = 16, 18, 20 were averaged for 
each g~ l , which should provide a good estimate of the asymptotic value, as can be seen in 
Fig. H The error was estimated based on the maximum deviation from this mean, among 
the three data points, taking into account the la error estimates that have been represented 
in the figure by error bars. Tabled shows that the large N expectation values (|4.4)l in the 
deformed theory are (up to statistical errors) the same as those of the undeformed theory. 

Incidentally, one might worry that due to the nontrivial moduli space of the undeformed 
theory, its partition function (for the purely bosonic theory) is not well-defined. Certainly 
this is the case for case for the Od reduction of d < 4 577(2) pure Yang-Mills; however, this 
is not the case — in spite of noncompact flat directions — for the Od reduction of d > 5 SU (2) 
pure Yang-Mills; see Eq. (24) of [12] . This is because, properly speaking, the classical moduli 
space is a set of measure zero in the field integration of the partition function. For d > 5, 
"entropic effects of the measure... overwhelm the possible divergences" ^2]. We are not aware 
of any exact result showing that the partition function of the CKKU quiver system obtained 
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5 10 15 20 25 



Figure 1: Trajectories of fixed lattice spacing, increasing volume, with data connected by 
lines to guide the eye. Each line is marked by the corresponding value of g -1 . A doubling of 
g^ 1 corresponds to a halving of the lattice spacing. 



g- 1 


(0*> ) 




2 


0.2469(94) 


0.494(19) 


4 


0.1249(27) 


0.500(11) 


8 


0.0618(22) 


0.494(18) 


16 


0.03064(75) 


0.490(12) 


32 


0.01586(64) 


0.507(21) 


1 (fi = 0) 




0.496(22) 



Table 1: Large N asymptotic values for the fi~ l = N trajectories. For comparison, the result 
for the undeformed (/i = 0) expectation value is shown in the bottom line. Estimated errors 
in the last 2 digits are shown in parentheses. 
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0.54 



0.52 
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0.46 




Figure 2: Undeformed theory results. Error bars differ in size due to differing numbers of 
samples. 

from the 6d — > Od SU(kN 2 ) pure Yang-Mills is divergent. 4 In our view it is unlikely, just from 
the result that 6d — > Od £77(2) pure Yang-Mills has a finite partition function. Furthermore, 
if there was a problem with the partition function of the undeformed theory, we would have 
expected an uncontrolled dispersion in our observables when we attempted to measure them 
in the undeformed theory. This is because flat directions not suppressed by entropic effects 
would allow for the Monte Carlo simulation to wander wildly throughout configuration space, 
leading to results that would not converge to reliable average values. That this did not occur 
is further evidence that the partition function is well-defined. Finally, the deformation Ssb 
can be regarded as a regulator of any possible divergence associated with the noncompact flat 
directions. That we obtain stable, identical results as the deformation is removed provides 
further evidence that our expectation values for the undeformed theory are reliable. 



5 Interpretation 

In our simulations of the purely bosonic system, we have observed that under the scaling 
([4.3)1 . the deformation becomes ineffective at changing the expectation value (|4.4jl . or equiv- 
alently (j4.1|) . away from the value that would be obtained in the undeformed theory. What 
has happened is that the flat directions that were lifted by the deformation are becoming flat 
all over again as iV — > oo. More precisely, the deformation is proportional to l/g^L 2 , and in 
the thermodynamic limit this quantity vanishes. Whereas the configurations of the moduli 
space of the undeformed theory that were lifted by Ssb cost energy at finite N, the number 
of such configurations is becoming vast due to the approximate flatness in those directions. 

4 The diagonal 17(1) from U{kN 2 ) = U{\) x SU{kN 2 ) decouples; its (divergent) contribution to the 
partition function trivially factors out. 
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For N ^ g 1 ' 2 , the entropy of these configurations wins out over the energy arguments that 
prefer (jSH|. 

One lifted region of the So moduli space where this is particularly clear is the branch 
()3.6|) . The action for such configurations is (setting c = 1) 



Integrating exp(— 5s) over all i:^ 3 we obtain 

N 2 



27 



cxp 



(-l+^ln^L 2 )) 



(5.2) 



Note that for large system size g2L ^> 1. Thus, the positive (entropic) term under the 
exponential wins out over the negative (energetic) term by a large margin. It would seem 
that the weight of these configurations increases exponentially as we increase N while holding 
g2L fixed; that is, in the continuum limit. 

Indeed the observed behavior summarized in Tabled is that in the deformed theory 

i Tr (x m x\Sj « |^ = ^g 2 a (5.3) 

Thus the simulations likewise indicate that configurations with x m = are dominating as we 
increase N while holding g 2 L fixed, which is nothing but the continuum limit. By symmetry, 
the same also holds for y m . 



6 Including fermions 

While the above results and arguments for the non-supersymmetric system are interesting 
in their own right, it is essential to include the fermions if we are to draw conclusions for the 
super symmetric system. 

Unfortunately, the present system suffers from a complex fermion determinant [2]. Thus 
the integration measure det Me~ Ss is not positive semi-definite and does not provide a 
satisfactory probability measure for Monte Carlo simulations. If we factor out the phase 
e ia _ j^j]^ anc j uge J ^g^. jVf| e -Sg instead, then we will generate the phase-quenched 

(p.q.) ensemble of boson configurations. Expectation values of an operator O in the full 
theory are formally related to those in the phase- quenched theory by the reweighting identity 

(O) = (e ia O) p J(en M . (6.1) 

In some lattice systems, the distribution of a in the phase-quenched ensemble is sharply 
peaked. The phase-quenched quantities in the ratio exist and can be measured reliably. For 
example, this is the case in the 4d U(1)l x U(1)r symmetric Yukawa model of [T^]. By 
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contrast, the results of 2 J for the CKKU system are not encouraging. The distribution of a 
(denoted in that work) in the phase-quenched ensemble was found to be essentially flat. 
This leads to approximate cancellations when we attempt to compute reweighted quantities 
contained in the numerator and denominator of (j6.1|) . More details will be given below. 
Here we merely emphasize that we do not expect the strategy implied by (j6.1J) to succeed, 
in practice, for the present system. Our chief result will thus be to quantify the degree to 
which it fails. 

As will be seen below in Sec. 16.21 det M vanishes for the Sq = configurations that 
were discussed above. Thus the configurations that have the most weight in the fully- 
quenched ensemble are configurations that have vanishing weight in the phase- quenched 
ensemble. (Further quantification of this statement is given in Sec. 16.31 ) For this reason we 
find it interesting to study the quantity (0 m °) p q , in Sec. 16.41 Certainly this bears a closer 
resemblance to the expection value in the full theory, in that the suppression due to very 
small | det M\ is taken into account for configurations in the neighborhood of the undeformed 
moduli space. In Sec. 16.51 we return to the matter of phase reweighting. 



6.1 Fermion action 

First we review the fermion action in the notation of |2|. The fermion action can be written 
in the form 



4 W 
9 



The elements of the fermion matrix are: 



pp 

S2,n 



(6.2) 



TUP 



(iVCJl.2 
« n )l,4 

(iVC n k4 



4,2 



(M, 



mn)2,2 



(M, 



mn)3,3 



+IJ.up v _j_ fHpv v 
l 'm,n^'n+i ' b xa,n^xii 

4.pVp—V _ 

_+PVp „ U -4- V-r- ... 

L m,n Jj n+j ' 1/ m,n+i J "ii5 



j-PP 1 * TfV 

t m,n-i a 'm> 
ppu 



(Ml)l,3 = 
« n ) 2l l = 

(MZ 



b m,n.yn+i 1/ m,n4-ji'n5 



.pup v 
m,n n+J 



1/ m,n x 'in 



2,4 



(/ m,nyn m,n— jomi 



(MZh,2 



.pvp —v 
t m,n x 'n 



,ppu —v 
(/ m,n / 'n) 



(MZ)4,1 



— +VP ryV . 



m,n+i J 'n' 
,ppv —v 



+P-vp —v 1 ./W — 1/ 



(Mi 



ran) 4,3 



1/ m,n' 6 n4-i4-J 



(6.3) 



7^ - 
t ooo 



;i 2 , a a ), Tr (T"T U T P ) = 2t pvp 



t aoo = 



«6 



icibc 



abc 



(6.4) 

This matrix has two ever-present fermion zeromodes. The first is associated with left 
multiplication (the transpose "T" merely indicates a row vector, for consistency with (|6.2)l ): 

(Ce2 P ,„,e 3 P n) A n ) r = (0,0,0,AO T V n (6.5) 
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The second is associated with right multiplication: 

Km , V«m , $U , XS.) = (0 , , , 0) V m (6.6) 

Because the matrix M is not hermitian, it is diagonalized by M — > D = U MV with U and 
V independent unitary matrices. When this is done, the diagonal matrix D always has 2 
zeros on the diagonal, one each from (|6.5|) and (|6.6() . 

The zeromode eigenvalues of the daughter theory can be factored out following the 
method used in El- We deform the fermion matrix appearing in (|6.2j) according to 

.\/ — .\/, - .\/ • el Nf (6.7) 

where Nf = 1QN 2 is the dimensionality of the fermion matrix and e C 1 is a deformation 
parameter that we will eventually take to zero. We factor out the zero mode through the 
definition 

M = lim M e , M e = e~ 2/N fM t det M = lim e" 2 det M t . (6.8) 

If this deformation is added to the action, it explicitly breaks the exact lattice super- 
symmetry and gauge invariance. This infrared regulator could be removed in the continuum 
limit, say, by taking ea <C iV _1 . Noting that L = Na is the physical size of the lattice, 
the equivalent requirement is that e <C L^ 1 be maintained as a — > 0, for fixed L. Thus 
in the thermodynamic limit {L — > oo), the deformation is removed. The parameter e is 
a soft infrared regulating mass, and is quite analogous to the soft mass \i introduced by 
CKKU P in their Eq. (1.2) to control the bosonic zeromode of the theory. In the same 
way that the deformation introduced with fi does not modify the final results of the renor- 
malization analysis of Section 3.4 of pQ, our e does not modify the result of the quantum 
continuum limit. The essence of the argument is that we have introduced a vertex that will 
be proportional to the dimensionless quantity g%ea 3 <C g%a 3 /L, where g 2 is the 2d coupling 
constant. Such contributions to the operator coefficients C, C in Eq. (3.29) of |1 vanish in 
the thermodynamic limit. Because the target theory is super-renormalizable, we are assured 
that the perturbative power counting arguments are reliable and the correct continuum limit 
is obtained. 

In 0, we studied the convergence of det M e — > detM - We found that the convergence 
was rapid and that a reliable estimate for det Mq can be obtained in this way. As a check, we 
computed the eigenvalues of the undeformed matrix M, using the math package Maple, for 
10 random boson configurations. We found that the product of nonzero eigenvalues agreed 
with det M in magnitude and phase to within at least 5 significant digits. 

In the phase-quenched ensemble, we take the probability measure defined by | det M t \e~ SB , 
with e = 10~ 6 . Our numerical studies indicate that this e is more than small enough to factor 
out the ever-present zeros while leaving the nonzero eigenvalues virtually unchanged. We 
have dropped the factor of e -2 that distiguishes between M e and M e . This is permissible 
because it cancels in correlation functions. In the simulation, this is reflected in the fact 
that under a shift in the boson configuration, ^lnldetMel = <51n|detM e |, since — 2 hie is 
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a constant. The value e = 10 -6 renders M e rather ill-conditioned, so we must use a very 
accurate inversion algorithm (see below). As an aside, we have found that approximation 
of M~ l by noisy estimators fails because the bicongujate gradient algorithm cannot handle 
the ill-conditioned matrix. This is the main obstacle that we have faced in our attempts to 
extend our analysis beyond N = Q. Note that for N = 6, the matrix M e is already 576 x 576 
with complex entries. 

6.2 Determinant on moduli space 

We have stated above that the effects of fermions must be taken into account in order to 
draw conclusions about observables of the supersymmetric system. One way to see this is 
to consider the fermion determinant in the case where we set the boson configuration to 
one of the points of the undeformed theory moduli space discussed above in Sec. 13.11 We 
now show that additional fermion zeromodes appear in these boson configurations. This 
seems reasonable if we regard the lattice action as the Hamiltonian of a classical statistical 
system. Zero action boson configurations correspond to bosons of zero energy. By the 
exact supersymmetry of the undeformed system, we expect that these zero energy bosonic 
states have fermionic partners, also of zero energy. These fermion zeromodes will cause the 
fermion matrix to be more singular. 5 Indeed, this heuristic argument seems borne out by 
the examples below. However, we were unable to constuct an explicit proof; presumably 
it would involve a study of eigenvalues of the transfer matrix corresponding to the lattice 
action. 

6.2.1 x m = y m = non-zeromode branch 

Here we can choose 



It is easily checked that M acting on this vector vanishes for any choice of the N 2 Grassmann 
variables A n . This gives N 2 fermion zeromodes associated with left multiplication. We also 
have N 2 fermion zeromodes associated with right multiplication: 



Note that the number of zeromodes, 2N 2 , is the same as the number of (complex) boson 
zeromodes found in Sec. 13.1.21 above. 

Two of these modes are the ever-present zeromodes discussed in Sec. 16.11 above. The 
other 2(N 2 — 1) are zeros special to this branch of moduli space. Strictly speaking, the set 
of configurations that satisfy x m = y m = is a set of measure zero in the lattice partition 
function. But taking a small neighborhood of radius 5 about this subspace — that is, small 
fluctuations about the classical configuration — we see that det M e oc e 2 8 2<yN ~ 1 ) in this region, 

5 These observations were pointed out to us by a referee of an earlier version of this report I 1 . 



(CCe3 P „,A^ n ) T = (0,0,0,A n O T 



(6.9) 



Mm » <,« > €m . Xm) = (0 , , ^ m 5°" , 0) 



(6.10) 
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where the powers of 5 come from the 2(N 2 — 1) fermion approximately- zero- modes that are 
present near the x m = y m = branch. Thus the inclusion of fermions dramatically changes 
the weight of configuration space in the neighborhood of x m = y m = branch in the 
integration measure; it has become totally negligible. Thus we do not expect the entropic 
effects discussed in Sec. [3] to persist in the super symmetric system. As a consequence, we 
expect to get a very different value for (<j) x '°). 

6.2.2 z m = non-zeromode branch 

First consider the case where u k ' 1+k ' 2 = 1. Then set 

(^»,Sn,eCn,AS) T = (0,O l ^,0) T V n (6.11) 

It is easy to check that this is a fermion zeromode on this branch, for any choice of £. There 
is also a corresponding fermion zeromode of right multiplication: 

W, m , $U , > Xm) = (0 , , , S^x) V m (6.12) 
Now consider the case where uj k '^ +k ' 2 = — 1. Then set 

(ei,n,e2, n ,e3,n,Kf = (o,o,s^(-r +n %or v n (6.13) 

It is only slightly more work to check that this is a fermion zeromode on this branch. There 
is also a corresponding fermion zeromode of right multiplication: 

W, m , , <m , Xm) = (0 , , , ^°(-) mi+m2 x) V m (6.14) 

We note that the suppression due to the fermion determinant is much smaller, relative to 
the x m = y m = branch, in the neighborhood of this branch: det M e cx e 2 S 2 . 

6.3 Determinant distribution in the fully-quenched ensemble 

Formally, the expectation value of an operator O in the phase-quenched ensemble is related 
to a ratio of expectation values in the fully-quenched (/.<?.) ensemble: 

(0) M = (| det M\0) f J (| det M\) f _ q _ (6.15) 

This is a type of reweighting identity. However, for this to work in practice, it is necessary 
that there be a significant overlap between the phase-quenched ensemble and the fully- 
quenched ensemble. We now show that this is certainly not the case in the present system, 
which is just to say that the fully-quenched expectation values have no bearing on the 
system with fermions. Our main point here is just to show that the inclusion of the effects 
of fermions in the updating process of the Monte Carlo simulation is crucial to the study of 
the supersymmetric system, a not very surprising result. We now quantify this statement 
by an analysis of the impracticality of ()6.15j) . 
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We accumulate N s points of data in the fully- quenched ensemble and then bin them 
according to In | det M\. From the number of samples n« in each bin i, we obtain an estimate 
W{ ~ rii/N s of the weight of each bin. We denote the average value of in each bin by 
(0m°)-i- We denote 

| det M\i = exp (central value of In | det M\ in bin i) (6.16) 

In principal, as the number of bins is taken to infinity and the resolution of each bin is taken 
to be infinitely fine, the approximations 

<|detM|^V « £ti*|detM| 4 <C°>i 

i 

(|detAf|) A . « ^Wi\detM\i (6.17) 

i 

become exact. In practice, however, we always have truncation errors from using a finite 
number of bins, and systematic errors from finite resolution. Furthermore, the statistics 
in all of the bins that contribute significantly to these sums must be large, so that error 
estimates for the quantites Wi and are reasonably small. 

In Fig.Elwe show Wi, an d Wi\ det M|j as a function of log 10 | det M\%. It can be seen 

that the overall weight Wi \ det M\i is sharply peaked in a region where Wi is very small. (On 
the plot, Wi is indistinguishable from zero in this region; numerically, it is (9(10 -4 ).) The 
estimate of (4>^ > ) i in these upper bins is very inaccurate, due to a small number of samples. 
In fact, even for 10000 samples, we were not able to see the right-hand side of the peak! Thus 
it is impossible to estimate the truncation error that results from having data for only a finite 
number of bins. For the present system, it is not possible to get a reliable estimate of either 
of the averages on the right-hand side of (|6.15|) . This is consistent with the results of Sec. 16.21 
the fully-quenched ensemble gives its largest weight to configurations that have practically no 
weight in the phase- quenched ensemble. Thus the overlap is essentially nil, rendering ()6.15|) 
quite useless. Another conclusion that we can draw from this is: expectation values obtained 
in the fully-quenched ensemble are not a good indicator, to any degree of approximation, of 
expectation values that would be obtained in the phase-quenched ensemble. 

6.4 Results for the phase-quenched ensemble 

In Fig.0](0^ ) is obtained in the phase- quenched ensemble; that is, the fermion determinant 
was taken into account in Metropolis updates, using the linear variation 

6 In | det M e \ = Re Tr [M e _1 <JM e ] (6.18) 

As usual, we worked at e = 10 -6 , so that the ever-present fermion zeromodes factor out. 
This renders M t ill-conditioned. To invert the matrix reliably, we have used triangular 
(LU) decomposition by Crout's method, which is very accurate. We have checked that the 
inversion is good to 1 part in 10 6 in our simulations. 
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Figure 3: The dashed line is u>j. The solid line connects data for Wi\ detM]j, marked by 
diamonds. Data marked by an "x" represent ((f>^) i - Where error bars are not shown for the 
data, only 1 or 2 points of data were available, so no estimate of error could be made. 
This plot is for 10000 samples on a 6 x 6 lattice with g^ 1 = 2. 
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Figure 4: Phase- quenched averages (lines to guide the eye). 



19 



1.00 



0.50 



0.05 



0.02 



- 






1 ' ' ' 


' 1 ' ' 


' ' 1 ' 


- 


- 






\ 






- 






















\ 














\ 

\ 

^ — — — — _ _ \ 














\ 

\ 








— 






\ 

\ 






— 




+ 


1/ga 2 


= 4, p.q. 










X 


1/ga 2 


= 4, f.q. 










© 


1/ga 2 


= 16, p.q. 










□ 


1/ga 2 


= 16, f.q. 

1 , , , 


, 1 , , 


, , 1 , 





Figure 5: A comparison of some of the averages obtained in phase-quenched (solid lines to 
guide the eye) versus fully- quenched (dashed lines to guide the eye) simulations. 

We remind the reader that larger values of g~ l = 1/ga 2 correspond to finer lattices, 
and thus extrapolate to the continuum limit. However, to hold L = Na fixed, successive 
trajectories must be compared as iV versus 2N, since g~ x is doubled. It can be seen from 
Fig. |3] that for N < 6 it is not possible to extrapolate to the large L behavior. The best we 
can say is the following. In Fig. El we provide a comparison between the phase-quenched and 
fully-quenched results. It can be seen that the effect of the fermion determinant is dramatic. 
For the finest lattice that we have available, g^ 1 = 16, the phase-quenched results are more 
than an order of magnitude larger when we look at iV = 6. Furthermore, the phase- quenched 
results are always of the same order of magnitude as the classical estimate (|3.2(J|) . 

6.5 Phase reweighting 

As has been reported earlier and mentioned above, the distribution of a = argdetM is 
essentially flat with respect to the phase-quenched ensemble. In Fig. HO we present results for 
a 2 x 2 lattice. It can be seen that it will be quite difficult to obtain a reliable estimate of 
(exp(ia)) M _. 

E.g., we find that for N = 2, g~ x = 2, with 1000 samples, the "noise-to-signal" ratio is 
ff <cM«) M ./(cosa) M . = a {siQa)pq /(sma) pq = 0(1) (6.19) 
For N = 4, g~ x = 2, with 1000 samples, we find 

°"<cosa) p . 8 ./(cosQ!) M = a (sina)p 9 /(sina) p (? = 0(10) (6.20) 
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Figure 6: Average frequency distribution for a = argdetM in the phase-quenched distri- 
bution for the 2x2 lattice. Error bars are determined by variance in bin counts for several 
blocks of data. 

On general grounds, these noise-to-signal ratios are expected to get exponentially worse as 
iV is increased. Similarly discouraging results are found for the other quantities that need to 
be estimated in order to make use of the reweighting identity (jfi.ljl . Since reweighting does 
not even work for small N, it looks like a futile approach by which to study the full theory. 

7 Conclusions 

In the fully-quenched system we have found that the classical estimate for the expectation 
value is not reliable. We have argued that this is due to entropic effects that dominate 

when the relative strength of the deformation is scaled to zero. Thus we conclude that the 
fully-quenched system does not yield the nonsupersymmetric gauge theory in the quantum 
continuum limit. However, we regard it as an interesting question whether or not a well- 
defined continuum theory is obtained if the relative strength of the deformation is not scaled 
to zero. We have not investigated this question as it is beyond the intended scope of this 
article. 

Next we have shown that the fermion determinant tends to vanish in the neighborhood 
of the configurations that were argued to yield the dominating entropic effects. It follows 
that the fully-quenched ensemble will not provide a reliable guide to the physics of the full, 
super symmetric theory. We have studied ((f) x '°) in the phase-quenched ensemble. We have 
found a dramatic increase in the expectation value that is obtained. We regard this as 
encouraging for the CKKU proposal. 
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However, to truly obtain data on the full lattice theory proposed by CKKU, it is necessary 
to take into account the complex phase of the fermion determinant. We have discussed the 
notion of reweighting by the phase of the determinant. However, we know from previous 
work that the distribution of this phase in the phase- quenched ensemble is so nearly flat that 
it is not practical to obtain reweighted quantities. As a practical matter, it seems impossible 
to study the continuum limit of the full theory by Monte Carlo methods, unless there is 
significant progress in the technique for the treatment of the complex phase of the fermion 
determinant. 

In this regard, we remark that the complex phase is not present in the target contin- 
uum theory; so, it may be possible to address the phase by blocking out the short distance 
modes of the lattice theory. However, this too is beyond the scope of this report. Alterna- 
tively, it may be possible to improve on the CKKU construction and eliminate the complex 
phase altogether. We are particularly interested in the constructions that do not involve a 
dynamical lattice spacing such as the one that appears here. Many attempts to formulate 
super symmetric field theories on the lattice may be found in the literature (see for example 

EH El EH E3 and references therein). In our opinion, the recent successes in non-gauge 
models are very encouraging |2~T| 12*2] . Of particular importance is the understanding of 
these models as "topological" or Q-exact, where Q is an exact supercharge of the lattice 
system I2H EH| Indeed, the CKKU undeformed action can be written in a Q-exact 
form. An exciting development has been Sugino's exploitation of this Q-exact idea to con- 
stuct lattice super- Yang-Mills with compact gauge fields and an ordinary (non-dynamical) 
lattice spacing j2H] . While the Sugino construction has its benefits, he notes (based on a re- 
mark by Y. Shamir) that these systems also suffer from a vacuum degeneracy problem that 
renders the classical continuum limit ambiguous. Consequently, Sugino has introduced a 
non-supersymmetric deformation in these models that lifts the unwanted vacua; he demands 
that the relative strength of this deformation be scaled to zero in the thermodynamic limit. 
Thus in some respects the Sugino construction has a feature that is similar to the models of 
CKKU. He argues that entropic effects do not destroy the vacuum selection imposed by his 
deformation. We are currently investigating various aspects of Sugino's proposal [2*7] . 
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A M = 4 moduli space 



Here we establish the well-known solution to ()3.3|) . One way to see this is as follows 
First we note that S reduced to the zero modes, which we write as S z , takes the form 



N 2 /1 

S z = ^Tr U[ X \ X ] + [y\y] + [ Z \ Z ]y 

+2[x, y] [y f , x f ] + 2[y, z] [z\ y f ] + 2[z, x] [x\ 



(a.i; 



Now note that the U(l) parts of x, y, z do not appear and can take any value. Thus we can 
restrict our attention the the SU(k) parts, which we choose to express in terms of Hermitian 
matrices a p , b p , p = 1,2, 3: 



x c T c 
y c T c 

z c T c 



(a£ + ibl)T c = a x + ib x 
{a c 2 + ib c 2 )T c = a 2 + ib 2 
{a c 3 + ib c 3 )T c = a 3 + ib 3 



(A.2) 



Substitution into (|A.1|) and a bit of algebra yields 



S z 



N 2 m 

g 2 



2 ^$^[ a f' 6 f] ) + 



[6p,a,])' 



iV 2 



— Tr ([op> S] 2 + h, Kf + 2[ap, &,] 2 ) 



(A.3) 



Using positivity arguments quite similar to those above, one finds that S z > and that 
S z = iff 



[op, a,] = [bp, b q ] = [a p , b q ] = 0, Wp, q 



(A.4) 



which is nothing other than Eq. (53) of j2Hl- Since the matrices are all hermitian and they all 
commute, it is obviously possible to choose a basis which simultaneously diagonalizes them. 
This basis will be related to the one used in (|A.2|) according to T c — > T' c = aT c a\ which is 
nothing other than the global gauge transformations (|3.4|) . 
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